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ABSTRACT 

I use the most recent cosmic microwave background (CMB) anisotropy 
measurements to constrain the leading cold dark matter models in the O m — £l\ 
plane. A narrow triangular region is preferred. This triangle is elongated in such a 
way that its intersection with even conservative versions of recent supernovae, cluster 
mass-to-light ratios and double radio source constraints, is small and provides the 
current best limits in the O m — £l\ plane: Oa = 0.62 ± 0.16 and Q m = 0.24 ± 0.10. 
This complementarity between CMB and other observations rules out fl\ = models 
at more than the 99% confidence level. 

Subject headings: cosmic microwave background — cosmology: observations 

1. INTRODUCTION 

The main goal of CMB measurements and the two new satellite missions MAP and Planck 
Surveyor is to determine a host of cosmological parameters at the unprecedented accuracy of a few 
percent (Jungmann et al. 1996, Zaldarriaga et aJ. 1997, Bond et al. 1997). As part of this goal it is 
important to keep track of what can already be determined from the CMB without conditioning 
on certain families of models or on certain values of parameters within these families. In this 
Letter, the analysis of CMB anisotropy measurements is expanded to include the most popular 
families of cold dark matter (CDM) models such as flat, flat- A and open models, as well as the less 
popular open— A models. Figure 1 presents an overview of this parameter space. Open-A models 
are considered here because they subsume all of the above models and thus provide a parameter 
space in which the most popular models can be directly compared. 

The popularity of non-zero S1a models has waxed and waned over the years; for excellent 
reviews see Felten <fc Isaacman (1986) and Carroll et al. (1992). Q\ was introduced by Einstein 
(1917) to solve the discrepancy between an apparently static universe and the dynamic cosmology 
of general relativity. Since this inauspicious beginning Q\ has been invoked many times and 
seems to be a surprisingly multi-purpose cure-all for theory-observation mismatches. Several 
recent papers (Turner 1991, Ostriker & Steinhardt 1995, Roos &: Harun-or-Rachid 1998), have 
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pointed out the effectiveness of 17a in resolving apparent conflicts between various observational 
constraints. Recently, S7a has been invoked to solve the discrepancy between globular cluster ages 
and the age of the Universe inferred from measurements of Hubble's constant. 

Recent supernovae results in £l\ = models yield Q m values so low that they are unphysical: 
Q m = —0.4 ± 0.5 (see column 2 of Table 1, error bars and limits in this Letter are 68% confidence 
levels unless stated otherwise). Not only are they unphysically low but the highest Q m values 
allowed by the error bars are in strong disagreement with the high values of Q m preferred by the 
CMB in these same models (Lineweaver 1998). In Lineweaver & Barbosa (1998b) we report a 
99.9% confidence lower limit of VL m > 0.31 (see Figure 1). This supernovae/CMB inconsistency is 
strong motivation to use CMB data to explore a larger parameter space which includes VLa- If the 
inconsistency is caused by the incorrect assumption that VLa = 0, then such an analysis will show 
it. The result of the analyis presented here is that U\ > can resolve this inconsistency. 

Testing large parameter spaces is important to minimize the model dependence of the results. 
For example, in Lineweaver & Barbosa (1998a) the CMB data favored h = 0.30^q'q7 (if VL m = 1, 
if VL\ = and if all the other assumptions made are valid). In Lineweaver & Barbosa (1998b), 
hereafter LB98b, we dropped the £l m = 1 assumption and still found low h values (h = 0.40 but 
with large error bars: 0.26 < h < 0.97) and Q m > 0.57. Thus the CMB data prefer high Q m values 
(if Qa = 0). These may be big if 's. The purpose of this paper is to make these if 's smaller by 
exploring a still larger region of parameter space. 

Other workers have used CMB observations to constrain cosmological parameters in CDM 
models ( e.g. Bond & Jaffe 1997, deBernardis et al. 1997, Ratra et al. 1997, Hancock et al. 1998, 
Lesgourgues et al. 1998, Bartlett et al. 1998, Webster et al. 1998, White 1998). The previous work 
most similar to this Letter is White (1998). White (1998) combined supernovae results with the 
Hancock et al. (1998) estimate of the position in I— space of the peak in the CMB power spectrum. 
In Section 5, I compare my results to White (1998) and other work. 



I use x 2 minimization to identify the best fit and Ax 2 = [1,4,9,16] to identify 
[68.3%, 95.4%, 99.7%, 99.9%] confidence regions around the minima. The x 2 computation 



where 9 = [Qa, ^m, h, ft^h 2 , n, Qiq] which are, respectively, the cosmological constant, the 
matter density, Hubble's constant, the baryonic density, the primordial power spectrum index 
and the primordial power spectrum amplitude at t = 10. f^A and Vt m are normalized to the 
critical density ($7 a = A./3H 2 , £2 m = p rn / Pcrit) and h is the dimensionless Hubble constant 
normalized at 100 km s^ 1 Mpc~ l . The following ranges for these parameters were used: 



2. 



Method 



is 




(1) 



-3- 



< ft A < 0.9, 0.1 < Q rn < 1.0, 0.15 <h< 1.00, 0.005 < Q b h 2 < 0.030, 0.58 < n < 1.42, 14.0 /iK 
< Qio < 23 /LtK with respective stepsizes 0.05, 0.05, 0.05, 0.005, 0.03, 0.5. Although I am exploring 
a 6 dimensional parameter space I limit the discussion in this Letter to the Cl m — tt\ plane. A 
X 2 value is calculated for all points in the ranges above, consistent with the condition Q b < Cl m . 
All 6 parameters are free to take on any value which minimizes the x 2 ■ See §2 of Lineweaver et 
al. (1997) and §2 of LB98b for more details of the method. Compared to LB98b the improvements 
here are: 1) inclusion of fl\ as an extra dimension of parameter space ( not just flat— A), 2) new 
and updated CMB data points, 3) a larger range and higher resolution in the f^/i 2 dimension, 4) 
the weak dependence of the helium fraction Y# e on Q^h 2 is included (Sarkar 1996, Hogan 1996). 
Previously we had just set Y# e = 0.24 for all values of tl^h?. 

The models in Eq. (1) were computed with CMBFAST (Seljak & Zaldarriaga 1996). The 
sum on i in Eq. (1) is over 35 independent CMB anisotropy measurements. These data points are 
given in Table 1 of LB98b with the following updates: two new points from Femema et al. (1998), 
updated values from Baker (1998) and Leitch (1998) to conform to the final published results and 
improved estimates of the error bars on the earlier MSAM results. The previous 5 MAX points 
have been combined into one point (Tanaka 1997) and I use the uncorrelated DMR points reported 
in Tegmark &; Hamilton (1997). I also now include the data point from Picirillo & Calisse (1993). 
In this analysis I use the Leitch (1998) recalibration of the Saskatoon measurements (Netterfield 
et al. 1997) which involves a 5% increase in amplitude with a 4% dispersion about the new central 
values (see Section 2.3 of LB98b). 

I have compared the results from LB98b with results from this updated data set in open 
models (Qa = 0). The difference is small. LB98b reported 0.26 < h < 0.97 and tt m > 0.53, 
(Q m > 0.31 at 99.9% confidence level) with central values of h = 0.40 and Q m = 0.85. The 
analogous result with the new data set is 0.30 < h < 0.98 and £l m > 0.57 (fi m > 0.34 at 99.9% 
confidence level) with central values h = 0.55 and Q m = 0.75. 

3. CMB Results 

The main result of the CMB analysis is shown in Fig. 2. The CMB data prefer the narrow 
hashed region (68.3% confidence level) enclosed by approximate 95.4%, 99.7% and 99.9% contours. 
The 'X' marks the best fit: (ft m , fi A ) = (0.45, 0.35). The minimum x 2 value is 22.1. With normally 
distributed errors, the probability of obtaining a x 2 this l° w or lower is 22.3%. 

Since we have not considered closed models, the 68.3% confidence region is cut off by the 
£l m + Qa = 1 limit. However the x 2 values along this cutoff are very close to the 68.3% confidence 
level. Thus when the entire square is explored one should expect the 68.3% confidence region to 
widen only for £l m ^ 0.5 and only by a narrow strip approximately parallel to the £l m + f^A = 1 
line. If one restricts the analysis to flat-A models, the best-fit is ^a = 0.1 with an upper limit 
ft A < 0.61. 
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This analysis also yields new constraints on the power spectral index and on the power 
spectrum normalization. At the x 2 minimum n = l.OG^g'JI and Qio = 18 ± 1 fiK. The best-fit 
Hubble's constant is h = 0.75 but with a large 68.4% confidence range 0.35 < h < 0.98. Thus h is 
not very usefully constrained by the CMB alone. The \ 2 minimum is obtained at the maximum 
value considered for the baryonic density: Q^h 2 = 0.030 with a 68.4% lower limit of 0.025. Thus 
the CMB prefers high f^/i 2 values. 

4. Combining Constraints 

Current observational constraints from supernovae, cluster mass-to-light ratios and double 
radio sources in the Q m — f^A plane are given in Table 1. To approximate the region of the 
Q m — plane favored by non-CMB observations, I form likelihoods from the limits in Table 1 
and from published contours (e.g. Riess et al. 1998, Fig. 6, Carlberg et al. 1998, Fig. 1, Daly et 
al. 1998, Fig. 1). I then form joint likelihoods C non - C MB = £sn x ^Clusters x £radio, where all 
terms are functions of Q m and The CMB results are combined with the non-CMB results in 
the same way: Ctot = C-cmb x C non -CMB- 

There are a variety of ways in which the limits can be selected and combined. My strategy 
is to be reasonably conservative by trying not to over-constrain parameter space. Practically 
this means using contours large enough to include possible systematic errors. For example, two 
independent supernovae groups are in the process of taking data and refining their analysis and 
calibration techniques. Table 1 lists their current limits (ref. 3 & 6) which are consistent. I 
combine each supernovae result separately with the CMB constraints and list the result in the 
row marked '+ CMB' directly under the supernovae result. However, the main result I quote in 
the abstract (last row of Table 1) comes from using the least-constraining of the two (ref. 3) in 
the conservative combination of non-CMB constraints, i.e., Ctot = Cqmb x C non _cM B(conservative) ■ 
Figure 3 shows each of the three terms in this equation. 

I apply the same strategy with the Carlberg et al. (1997, 1998) cluster mass-to-light ratios. 
They report 30% errors in their Q m result but also cite a "worst case" of a 73% error if all 
the systematic errors conspire and add linearly. In Table 1, I give the result of combining the 
30% and 73% versions separately with the CMB constraints. However, I use the 73% error 
in the conservative combination of non-CMB constraints. Thus the ^a = 0.62 ± 0.16 and 
VL m = 0.24 ± 0.10 quoted in the abstract are conservative in the sense that the error bars from the 
SNIa and cluster mass-to-light ratios are "worst case" error bars. A summary of the systematic 
error analysis of the cluster mass-to- light ratios result is given in Section 9 of Carlberg et al. (1997) 
and of the supernovae results in Section 5 of Riess et al. (1998). 

The conservative result we quote is robust in the sense that when any one of these non-CMB 
constraints is combined with the CMB constraints, similar results are obtained: S7a = is more 
than the 95.4% confidence level away from the best-fit. Systematic errors may compromise one or 
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the other of the observations but are less likely to bias all of the observations in the same way. 

Those confident in the Riess et al. (1998) and the Carlberg et al. (1997, 1998) results 
should quote the extremely tight limits labelled 'optimistic' in the penultimate row of Table 1: 
VL m = 0.18 ± 0.04 and ft A = O^ll^g. In this small region of the VL m - Q A plane the CMB data 
prefer h = 0.80 ± 0.10, Q w = 18 ± 0.5 /jK, n = 1.0 ± 0.1 and Q b h 2 ~ 0.025. 

It would be useful to add constraints on £l\ from lensing. Although Kochanek (1996) and 
Falco, Kochanek & Muhoz (1998) report < 0.7, the lensing estimates from Chiba & Yoshii 

(1997) (J7a ~ 0.8) and Fort et al. (1997) (Q\ ~ 0.6) are in agreement with the result found here. 
Thus lensing limits on Qa still seem too uncertain to add much to the analysis. However, when I 
fold Qa < 0.7 into the analysis the result for comes down less than 0.5 a. 

5. Comparison with Previous Results 

White (1998) combined his supernovae analysis with CMB results and pointed out the 
important complementarity of the two. He used the Hancock et al. (1998) estimate of the position 
in £— space of the peak in the CMB power spectrum based on £l (= f2 m + ^a) — dependent (but 
not fi a - dependent) phenomenological models introduced by Scott, Silk & White (1995). The 
Hancock et al. (1998) result is: £ pea k = 263lg4 9 . This should be compared to the LB98b result of 
ipeak = 260^20- The tighter limits are presumably due to the more precise parameter dependencies 
of the power spectrum models and the more recent data set. In Figure 3 of White (1998), the la 
contours from supernovae and CMB overlap in a region consistent with the results reported here. 

Another result of the Hancock et al. (1998) analysis is Q m + VL^ = 0.7*qJ. This is consistent 
with the CMB results presented here. A rough approximation of the elongated triangle in Figure 
2 is a strip parallel to the Q m + S7a = 1 line approximated by VL m + f^A = 0.8 ± 0.3. 

By combining CMB and IRAS power spectral constraints, Webster et al. (1998) obtain 
£l m = 0.32 ± 0.08 in flat-A models with n = 1. This is in very good agreement with the values 
obtained here from the combination of CMB, supernovae, cluster mass-to-light ratios and double 
radio sources. 

The result presented in this Letter is consistent with a large and respectable subset of 
observational constraints (Turner 1991, Ostriker & Steinhardt 1995, Roos & Harun-or-Rachid 
1998). However models in this region of the Q m — f^A plane appear to have problems fitting the 
shape of the large-scale power spectrum measured from the APM survey around k = 0.1 h Mpc^ 1 
(Peacock 1998) and are in disagreement with constraints on VL m from the POTENT analysis of 
the local velocity field (Dekel et al. 1997). 

We have assumed Gaussian adiabatic fluctuations. However Ferreira et al. (1998) have 
analyzed the DMR four-year maps and found tentative evidence for non-Gaussianity. Peebles 

(1998) has presented the case for isocurvature rather than adiabatic initial conditions. If either 
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Table 1: Non-CMB and Non-CMB + CMB Constraints on 17 A and fl 



Reference) Method 


n A = o 


fi m + Ha = 1 


n m + n A < i 








n A 




n A 


1) SNIa 

2) SNIa 

3) SNIa 

3) SNIa + CMB 


U - 8S -0.60 
-0.2 ±0.4 

-0.4 ±0.51" 


94+ u - 34 

u - a -0.28 

0.6 ±0.2 
0.27 ±0.3 6 


0-06^0.34 

0.4 + 0.2 
0.73 ±0.3 


Q OO + 0.15 


^2+°' 25 


4) SNIa 

5) SNIa 

6) SNIa 

6) SNIa + CMB 


-0.1 ±0.5 

-0 2+ 1 - 
u - -0.8 

-0.35 ±0.18 


0.65 ± 0.3 C 
4+0-5 

24 +0.17d 


0.35 ±0.3 
0.761°;!? 


26+ 09 


63+ 14 


7) Cluster M/L 

7) Cluster M/L + CMB 

7) Cluster M/L 

7) Cluster M/L + CMB 


0.19 ± 0.14 e 
0.19 ± 0.06 / 






0.24 ±0.10 

17+ - 04 
V - L 1 -0.03 


62+° 17 
7 o+0.07 


8) Double Radio Sources 

8) Double Radio Sources + CMB 


-0 l +,) ' 5 


2+ u ij 

u - -0.2 




33+ 16 


0.52 ± 0.25 


Non-CMB(optimistic) 9 
Non-CMB(conscrvative)' 1 
Non-CMB (optimistic) + CMB» 
Non-CMB (conservative) + CMB h 








0.16 + 0.05 
0.15 ±0.13 
0.18 ±0.04 
0.24 ±0.10 


0.57 ±0.25 

°- 46 -0.41 

71+ - 07 

u -' i -0.08 
0.62 ±0.16 



a —0.4 ± 0.1 zb 0.5 (statistical and systematic errors respectively) I have added the statistical and systematic errors in quadrature 
b 0.27 ± 0.06 ± 0.3 (statistical and systematic errors respectively) I have added the statistical and systematic errors in quadrature 
c as quoted in Ricss ct aJ. 1998 

d Ricss ct aJ. (1998), Figure 6 ('MLCS method 1 + 'snapshot method') using either the solid or dotted contours whichever is larger (corresponding 

respectively to the analysis with and without SN1997ck). 

e "worst case" result with 73% errors (Carlberg et ai. 1997, p 473) 

^ 30% error cited as the main result 

9 optimistic combined constraints using SNIa results from ref. 6) (rather than rcf. 3) and using the 30% error bars (rather than the 73% error 
bars) from Carlberg ct aJ. (1997) 

h conservative combined constraints using SNIa results from rcf. 3) (rather than rcf. 6) and using the 73% error bars (rather than the 30% error 
bars) from Carlberg ct aJ. (1997) 

1) Pcrlmuttcr ct ai. (1997), 2) Pcrlmuttcr ct ai. (1998), 3) Pcrlmuttcr private communication (1998), 4) Garnavich ct aJ. (1998), 5) Schmidt ct 
a/. (1998), 6) Ricss ct ai. (1998), 7) Carlberg ct aJ. (1997) and Carlberg ct ai. (1998), 8) Daly ct ai. (1998) 
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non-Gaussian processes or isocurvature initial conditions play significant roles in CMB anisotropy 
formation, then the CMB results presented here are significantly compromised. 

6. Summary 

The results presented here are largely observational but are model dependent. In a series of 
papers (Lineweaver et al. 1997, Lineweaver & Barbosa 1998a, 1998b), and now in this work, we 
have looked at increasingly larger regions of parameter space. Each time the \ 2 minimum has 
been found to lie within the new region. This might be taken as a sign of caution not to take the 
currently favored region too seriously. On the other hand, our choice of new parameter space to 
explore has been guided by independent observational results. 

I have used the most recent CMB data to constrain the leading CDM models in the Q m — Qa 
plane. A narrow triangular region is preferred. This triangle is elongated in such a way that its 
intersection with even conservative versions of other constraints is small and provides the current 
best limits in the plane: f^A = 0.62 ± 0.16 and 0, m = 0.24 ± 0.10. This complementarity between 
CMB and other observations rules out f^A = models at more than the 99% confidence level. 

Until recently observations could not discriminate between a zero and a non-zero cosmological 
constant. However a wide range of observations have indicated that O m < 1 and the most recent 
observations appear to favour Qa > 0. The addition of the CMB constraints presented here to 
these other cosmological observations strengthens this conclusion substantially. 

I gratefully acknowledge discussions with Saul Perlmutter and Brian Schmidt about the 
supernovae data sets and I acknowledge Uros Seljak and Matias Zaldarriaga for providing the 
Boltzmann code. I am supported by a Vice-Chancellor's fellowship at the University of New South 
Wales, Sydney, Australia. 
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Fig. 1. — Overview of cold dark matter (CDM) models in the Q m — Q\ plane. This parameter 
space includes the popular cosmological candidates: flat, flat-A and open models as well as the less 
popular open— A models. Closed models (upper right triangle) are not considered here. The white 
contours indicate the region favored by the joint likelihood of constraints from recent supernovae, 
cluster mass-todight ratios, and double radio sources. They are approximate 38.3%, 68.3% and 
95.4% confidence regions (see Table 1 and Section 4). In open models ($7 a = 0) there is a significant 
inconsistency between CMB results and these other observations. Previously reported constraints 
from the CMB in open models are indicated just below the £l m axis: £l m = 0.85 with a 68.3% lower 
limit of 0.53 and a 99.9% lower limit of 0.31 (Lineweaver & Barbosa 1998b). 
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Fig. 2. — The current CMB data prefer the narrow hashed triangular region (68.3% confidence 
level) enclosed by approximate 95.4%, 99.7% and 99.9% contours. The 'X' marks the best fit: 
({l m , Qa) = (0.45, 0.35). This region is bounded by the limits: S7a < 0.77 and Q. m > 0.15. Assuming 
= 0, yields the lower limit O m > 0.57. Restricting consideration to flat-A models, the best-fit 
is fl\ = 0.1 with an upper limit Q\ < 0.61. 
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Fig. 3. — The thick dark lines are the approximate 68.3%, 95.4% and 99.7% contours from the 
joint likelihood of the CMB and non-CMB constraints. The 68.3% and 95.4% confidence levels 
from non-CMB observations are in white (same as in Figure 1). The 68.3% and 95.4% confidence 
levels from CMB observations are the thin black lines (same as in Figure 2 but partially obscured 
here). The best fit is : f2 m = 0.24 ± 0.10 and Q\ = 0.62 ± 0.16. = models are excluded at 
more than the 99.7% confidence level. 



